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intrinsic autoregressions have not been known. Various asymptotic expan- 
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is derived. The expression involves Appell's hypergeometric function of the 
fourth kind. Various particular cases of the expression are also derived. 
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1. Introduction 

Let {Xu,v : w, w G Z} be a homogeneous first-order intrinsic autoregression on 
the two-dimensional rectangular lattice [TT], [3] with generalized spectral 
density function 

/(w,?7) = K (1 — 2a cos w — 2bcosi])~^ 
for w e (— TT, tt), 77 G (— TT, tt) and the conditional expectation structure 

) + b {Xu,v-1 + Xu,v-\-l) (1) 

where a>0,b>0, a + b~^ and 

var {Xu,v I ■■ ■) = K. 

We can assume without loss of generality that k — 1 and that {X„,„} is Gaus- 
sian, [4]. It follows that the difference Xu,v — X^+s.v+t has a well defined distri- 
bution with zero mean and lag (s,i) variogram 



1 r r 1 - cos(sx) cos(iy) 
""^ Jo Jo 1 ~ 2a cos X — 2b cos y 
1 



u,,ia,b):=- I I , ^:Z^ d^dy- (2) 
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The computation of ([2]), both analytically and numerically, has been a subject 
of considerable interest. The symmetric case a — b = 1/4 was considered by 
McCrea and Whipple pSj , Spitzer [T^ and Besag and Kooperberg [3] . Besag and 
Mondal [4] derived explicit expressions for Vg^s and h'g^o for the asymmetric case 
with the latter expressed as a finite sum of incomplete beta functions. Duffin 
and Shaffer [9l Theorem 4] and Besag and Mondal f^*, Theorem 2] provided 
asymptotic expansions for v^^t in terms of r = ^/4bs^ + 4at^ for the symmetric 
and asymmetric cases. Their approach was to exploit the recurrence equation 
with respect to integer time variables s and t ^ Eq. (4)]. See also [2] and [TT| . 

The aim of this note is to derive an explicit expression for ([2|) for general 
s, t, a and b. Our approach is quite different — we attack the problem varying 
a and b to find the explicit formula for Vst{cL, b). 

The expression given in Section 2 involves Appell's hypergeometric function 
of the fourth kind [lOl page 1008] defined by 

^4«.ft7,V;..,]^ffMffi; (3, 

for \/\x\ + i/jyl < 1, where {w)i := w{w + 1) • • • (w + ^ — 1) denotes Pochham- 
mer symbol with (w)o = 1. Various particular cases of the general expression, 
involving simpler functions, are derived in Section 3. 

A transformation formula between and the Appell's hypergeometric func- 
tion F2 defined by 

J? \ R R' ' 1 (a)j-+fc(/3)j-(/30fc 11,11/1 
F2[a,(3,(3-nn;x,v\ :=2^2^ , , \x\ + \y\<l 

j=0 k=0 J- 

is [13 §5.4] 

F4[a,(a + l)/2;7 + i,7' + i; x^y^] 

r 2x 2v 

= (l + a; + j/)-"^^2 a,7,7';27,27'; ^ 



x + y + 1 X -\- y + 1 

This formula has been proved earlier (in equivalent form) in ^ p. 11, Eq. (3.1)], 
see [m Eq. (175), §9.4.] too. In-built numerical routines for the computation 
of F2 are available, see [H §3.1.3] and especially [7| §2]. Reduction procedures 
for ^4 to lower order analytical expressions have been developed by Niukkanen 
[Tl] , Paramonova and Niukkanen [13] and references therein. 



2. Main result 



Theorem 1 expresses ^ in terms of Appell's hypergeometric function of the 
fourth kind defined in ([3]). It applies for \a\ + \b\ < i. The case \a\ + \b\ = ^ is 
considered by Theorem 2. 
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Theorem 1. For all s,t G No, \a\ + |6| < ^ we have 



s 



X Fil^is + t + 1), ^{s + t) + 1: s + l,t + l:Aa^ , ib^]. (4) 
Proof. By equation 3.915(2) in [10], we can write 



fir /'TT /'OC 

-(1 — 2a cos x — 2b cos y)z 



dxdydz 



1 f f 

Vst{a,h) = — III {l ~ cos{sx) cos{ty)) ' 
Jo Jo Jo 

= ^ e""^(jo(-2aa;i)Jo(-26xi)-i"+*J4-2axi)Jt(-2&a;i)) dx (5) 

where i = \/— T and Jg (a;) denotes the Bessel function of the first kind of order 
s. Setting 

/•OO 

Ist{a,b) ■.^1"+* / e'^Js{-2axi)Jt{-2bxi)dx 
Jo 

we can rewrite ([5]) as 

Vstia^b) ^lQo{a,b) ~Ist{a,b) . (6) 
Consider now for instance the Laplace-transform formula 3.12.15(20) in jTB] . 



VIZ 



e-PH^ Jf,{At)J^{Bt) dt 



2M+^pA+M+<>+i r(^ + i)r(i/ + 1) 



X Fa 



i(A + ^ + i/ + l),i(A + At + /^ + 2); -, 



(7) 



when ^{\ + n + v} > -1, 3?{p} > \^{A}\ + Setting p = 1, A = 0,A = 

— 2ai, B = — 26i in ([7]) one deduces 

Ist(a, ^) = ^ ^) a'^* -P^4 [5(5 + i + 1), 5(5+) + 1; s + 1, t + 1; 40^, 46^] . (8) 

The result of the theorem follows from ^ and ([5]). □ 

Remark 1. The explicit expression in ([4]) has some applicability. Although one 
assumes that \a\ + \b\ = ^ in ^ the case \a\ + |6| < ^ has been of interest. 
For instance, see Example 3.3] quoting a result by Kempton & Howes on 
a four-neighbour auto-normal scheme with the maximum likelihood estimates 
a = 0.4848 and 6 = 0.0132, i.e. S + 6 = 0.4980 < \. 

Theorem 2. For all a E (O, 5) we h 



ave 



Vst{a,^-a)= lim iyst{aVT~e,{^- a) VT~9) . (9) 

— >0-\- 
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Proof. Consider the Laplace transform ([7]) for p = (1 — 9)^^^^ and A = 0. The 
resulting Appell's function F4 obviously converges on the edge a + b ~ ^ of the 
convergence region. Considering that Vgt is a difference of two terms, the 
result follows by Abel's summation method. □ 

Remark 2. The limit in Q gives an approximation formula, viz 

\ 2 , 



i.,t(a,i-a)«F4[i,l; 1,1; Aa\l ^ e),A{l - a) {1^9 

+ t 



s J (l-6»)-('*+*)/2 



F4 



s + l,t+l;Aa^{l- 9), 4(i-a)^(l- 



(10) 



where the quality of the approximation is controlled by suitably chosen small 
9>0. 



3. The symmetric case a — b — 1/4 



Here, we present technical details to calculate i/g 
formula [STJ §9.4, Eq. (149)], we can transform 

F4[a,(3;-f,j'; x,x] = 



4) 



Using Burchnall- 



a, (3, i(7 + 7'), 5(7 + 7'-!) 
7, 7', 7 + 7' - 1 



4a; 



The asymptotics of generalized hypergeometric series p+iFp,p > 3 have been 
studied by Biihring and Srivastava [B], Saigo and Srivastava [191 and by A.K. 
Srivastava i20|. Since we are faced with 



4F3 



i(s + < + l), i(s + t) + l, i(s + f) + l, i(s + t + l) 
s + 1, t+1, s + t + 1 



(11) 



we have a hypergeometric term 



4F3 



ai, a2, as, 04 



' 1 



that is zero-balanced, i.e. when ai + Q!2 + as + 04 = /3i + /32 + /^a- To account 
for the asymptotics of this ^F^l- • • ; 1 — 6*] as — > 0+, consider equation (4.2) 
in for p = 3: 



ai, a2, as, 



Pi, P2, 133 

valid ior 9 ^ 0+ . Here 



{L + B-ln9} +0{9) + 0{9\n9), (12) 



r(Q;i)r(Q;2)r(a3)r(a4) 
r(/?i)r(/32)r(/33) 



(13) 
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and 

where 7 = — ^'(1) denotes the Euler-Maschcroni constant and ^{x) — ^ Inr(x) 
denotes the digamma function. The B in p2p can be represented by various 
forms, see (4.6), (4.7), and (4.12) in 6 . An example is [6, Eq. (4.7)]: 



fe=i 

X 3F2 



(fe + /?! - Q4 - Q3)fc(/33 + 02- aj- a3)k 
k{ai)k{a2)k 

- as, (33 - q;4, -k 
P3 + Pi - a4 - a3, /33 + /32 - (34 - as 



1 



(14) 



An advantage of this formula is that it is a single infinite series of hypergeometric 
terms and the series for each hypergeometric term is finite because of the — fe in 
the numerator. Combining (fTTjl . (fT3|) and (fT4|) . we can write Lst ■= L, Bst ■= B 
and Tst ■= T as 

4^-(s+i)jt+i: 



2Jk 



k=l 



k{s + t+l)2k 



■ 3F2 



\{s + t), i(s + i + l), -k 



1 



(15) 



and 



S + t\ TT 

S 



Tst=T = 

respectively. Note for instance that 

B,, = *(.s + l)-*(s + i) seNo: 
see the Appendix. Finally, routine calculations show that: 

'l^ 4 ' 4 



(16) 



In 4 + Bn 



In 61 



s + 1\ (1 
s 



45+ 



'— { L,t + Bst - In [ + 0{9) + 0(6 In 9) 

^ St ^ 



= i{lnl6 + Boo + 27 + *(^(s + t + l)) + *(i(s + i) + l) -B,t} 
+ 0(61) +0(61 In 6*) 9^0+. (17) 

for all (s, t) G N§. 

Theorem 3. For all s,t e Nq we have 



Vst 



(18) 



Proof. In ([TT]) assume that s + t is even. So, applying properties of digamma 
function ^ to expression ([T7]), we conclude 

*(i(. + t+l))+vI/(i(,s + t) + l)=-27-ln4+ ^ +^ 



fc=0 



fc=i 



For odd s + t repeat this procedure. 



□ 
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Corollary 3.1. For all s G No, 



'1 1^ = 4 ^ 1 

k=0 



(19) 



Proof. Consider HH) for s = By (HH) it follows 



.(ii)^^{ln4 + 7 + *(.s + i)} = ^|ln4 + 7 + *(i)+£^| 



so the result. 



□ 



Remark 3. The equation (|19p is well-known in the literature, see [3] Example 
3.2], where the authors refer back to \18\ p. 148]. 
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Appendix 



Here we prove (IT6l) . Using the transformation given by Biihring O Eq. (4.1)], 
one can express Bst in pS]) as 



3F2 



i(s-)-t), i(s-t), 



i i(s-i+l)-A: 



1 



If s = t then the hypergeometric term reduces to 1, so 



]^,k{s + l). 



1 00 ^,^^^k-^ 



E 



E 



dx 



2(5 + 1)^0 tl (s + 2)fc (fc-1)! 



da; 



2(s 



TT)/ ^""4 J;2 ; 



dx 



r 3 



2(s + 1) 



3^2 



1, 1 

2, 2 



= *(s + l) -*(,+ !), 

where the last two steps follow by equation 1.512(5) in [TU] and equation 
3.13 (42) in [12], respectively. The proof is complete. 
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